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Final  Report 

Grant  Ho.  AFOSR-84-03 1 1 

"Development  of  Symbolic  Computation  Methods  for  Nonlinear 
Dynamics" 

Principal  Investigator:  R.H.  Rand 

Dept.  Theoretical  &  Applied  Mechanics 
\  Cornell  University,  Ithaca  NY  14853 

\  607-255-7145 

This  grant  was  used  to  purchase  the  following  equipment: 
SYMBOLICS  3670  Mainframe  Computer  (including  6  Kbytes  of  memory 
and  a  470  Mbyte  hard  disk) 

Laser  Graphics  Printer 
Ethernet  board 

and  the  following  associated  software: 

MAC3YMA 

FORTRAN 

PASCAL 

The  equipment  has  been  used  to  facilitate  computation  on  the 
fallowing  projects: 

1.  Normal  Forms 


Under  the  supervision  of  principal  investigator  R.H.  Rand, 
software  has  been  written  in  MAC3YMA  which  automatically 
performs  normal  form  computations  for  systems  of  nonlinear,  ' 
nonautonomous  differential  equations.  Together  with 
postdoctoral  research  associate  U.L.  Keith,  we  have  produced  a 
package  which  permits  the  user  to  perform  Taylor  expanded 
near-identity  transformations  with  unevaluated  coefficients  or.  a 
system  of  autonomous  CDE’s  (valid  to  terms  of  arbitary  order), 
and  then  to  choose  the  transformation  coefficients  so  that  the 
resulting  system  is  in  normal  form.  A  summary  of  this  work, 
including  the  program  listing  has  been  published  [1]. 

This  work  has  been  applied  to  the  the  nonlinear  parametric 
stiffness  control  of  flexible  systems  by  Professors  Moon  and 
Rand  [2],  and  to  the  dynamics  of  coupled  van  der  Pol  oscillators 
by  graduate  student  T.  Chakraborty. 

This  vxork  is  currently  being  extended  to  cover  a  class  of 
nonautonomous  periodic  systems  by  graduate  student  J.  Goon.  In 
this  case  the  computer  algebra  code  performs  both  Taylor  and 
Fourier  expansions. 

2.  Determinacv  of  Degenerate  Equilibria 

We  have  also  used  MACSYKA  on  the  3670  to  study  the  problem 
of  the  determinacy  of  degenerate  equilibria  with  linear  part 

(1 )  x'  =  y,  y*  =  0. 


Takens  [3]  has  shown  that  any  nonlinear  autonomous  system  with 
such  a  linear  part  can  be  put  in  the  normal  form: 
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is  topologically  equivalent  to  that  of  the  truncated  system: 
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(4)  x 1  =  y ,  y’  =  a2  x  . 


Ue  extended  his  result  in  the  case  that 
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E.g.  we  shoved  that  if  a?  =  0  but  a^  ^  0,  then  the  flow  (2)-(3) 
is  equivalent  to  the  truncated  system: 
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The  procedure  involved  transforming  from  cartesian  to  polar 
coordinates  in  the  neighborhood  of  the  singular  point  and  Taylor 
expanding  the  polar  flow  (called  blowing  up  the  singularity), 
then  testing  for  determinacy.  If  the  resulting  system  is  not 
structurally  stable,  then  we  blow  it  up  again.  In  some  cases 
this  required  9  consecutive  transformations  and  Taylor 
expansions,  a  computation  which  could  not  have  been  done  without 
a  KACSYMA  program.  This  work  will  soon  appear  in  print  [4]. 

3.  Hilbert’s  16*b  Problem 

Professors  Guckenhe  imer  and  Rand  have  been  using  HACSYIIA  to 

approach  Hilbert's  16th  problem  from  a  bifurcation  point  of 
view.  The  problem  is  to  find  the  maximum  number  of  limit  cycles 
exhibited  by  a  system  of  2  autonomous  OBE’s  with  quadratic  right 
hand  sides: 

2  2 

(7)  x'=ax  +by  +  cxy  +  dx  +  ey  +  f 

2  2 

(8)  y'  =  A  x  +  By  C  x  y  +  D  x  +  E  y  +  F 

This  problem  has  a  long  history  including  a  paper  by 
Petrovsky  and  Landis  [5]  which  supposedly  showed  that  the 
maximum  number  of  limit  cycles  was  3,  but  was  then  claimed  to 


have  errors  in  it  by  Moser  [6].  The  question  was  illuminated  by 
the  recent  example  of  Songling  [ 7 ] »  which  exhibits  4  limit 
cycles. 

We  have  been  using  bifurcation  theory  to  investigate  this 
problem  in  the-  following  way:  We  look  for  the  standard  type  of 
bifucations  which  give  rise  to  limit  cycles  (Hopfs  and  saddle 
connections),  but  require  the  bifurcations  to  be  as  singular  as 
possible.  E.g.'  in  the  case  of  the  usual  Hoof  bifurcation,  a 
single  condition  generates  a  single  limit  cycle.  However,  we 
have  shown  that  requiring  2  additional  conditions  will  give  a 
bifucation  involving  the  birth  of  3  concentric  limit  cycles. 

This  computation  pushes  the  3670  to  its  limits,  and  certainly 
could  not  have  done  on  a  smaller  machine.  Some  algebraic 
expressions  which  result  are  15  Kbytes  in  length. 

4.  High  Precision  Numerical  Work  on  MACSYMA 

Professors  Parlange  and  Rand  have  worked  on  a  problem  in 
scattering  theory  [8]  which  involved  a  slowly  convergent  series. 

50 

The  leading  terms  were  of  order  10  although  the  series  summed 
to  a  value- of  order  unity.  In  order  to  keep  track  of  the 
differences  between  successive  terms,  we  used  MACSYMA’s  ability 
to  work  in  arbitrary  precision.  We  summed  170  terms,  keeping 
100  digits  for  each  term.  Note  that  even  quadruple  precision  in 
FORTRAN  would  not  have  been  sufficient  here. 

5.  Liapunov-Schmidt  Method  in  Dynamics  Problems 

Professor  Rand  and  Dr.  D.  Armbruster  have  worked  on 
formalizing  the  Liapunov-Schmidt  method  to  problems  involving 
Hopf  and  other  bifurcations  in  systems  of  nonlinear  partial 
differential  equations.  This  method,  based  on  the  Fredholm 
alternative  theorem,  offers  an  alternative  to  traditional 
perturbation  methods  such  as  Linstedt 1 s  method  (cf.  [9])  or  the 
method  of  averaging. 
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